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This report is a compilation of research designed to evaluate the effec- 
tiveness of the Martinec-Rango model in laige watershed applications. In 
addition, research was conducted into the application of the Martinec-Rango 
model to real time forecasting situations. The model as described by Rango 
and Martinec (1979) may be stated in the following form: 


Qn+l ■ c n ( I*An< T A> Sta + PAn) (2) 

' 12 

86400 

+ l*Bn( T B) s Bn + p BnJ A B* (5280)^ 

12 

86400 

♦ t a Cn< T C> s Cn + p Cnl A C* (5280) 2 \ (l - k^) + Q^nfi 

86400 ' 


where 


0 

Cn 

*n 

T A 

®An 

p n 

A 

A,B,C 

(5280) 2 

12 

86400 “ 


is the average dally discharge [ft-*s~*] 
is the runoff coefficient 

is the degree-day factor [in* °F“^ * d“*J 
is the calculated number of degre-days in Zone A 
is the snow coverage (100Z - 1.0) 
is the precipitation contributing to runoff [in] 

is the area [mi 2 ] 

is the recession coefficient 

is an index referring to the sequence of days 
as subscripts refer to the three elevation zones 

converts in • mi 2 per day to ft*^s“* 


Martinec and Rango (1979) documented results applying a relatively simple 
deterministic model for making snowmelt runoff using satellite derived snow 
cover as a critical input parameter. The so-called Martinec-Rango model 
appeared to be promising as a potential forecast tool but it was believed 
that further testing on other watersheds was desirable to define the 
applicability of the model to various size and shape watersheds. In 
response to this need the USDA Soil Conservation Service and Resource 
Consultants, Inc., entered into a cooperative agreement with the National 
Aeronautics and Space Administration (NASA) to conduct further trials on 
the model's performance and make modifications as necessary. Two watersheds 
in southwestern Colorado were chosen as test watersheds. These were the 
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South Fork of the Rio Grande River and the Conejoa River. The results of 
these investigations have been reported by Shafer et al. (1981a) and by 
Jones et al. (1981). The results of the simulation shoved the model 
capable of explaining an average of 89 percent of the observed streamflov 
variation on the South Fork of the Rio Grande and 87 percent on the Conejos 
and for the period 1973-79. Streamflov simulations vere run for the April 
1-September 30 interval. 

As an outgrovth of these previous investigations, further testing of the 
model vas conducted on a large vatershed (greater than 1,000 square miles). 
In addition, various techniques vere studied to develop a real time fore- 
casting version of the Martinec-Rango model. This report summarizes the 
results of these further studies. 

These additional studies on the Martinec-Rango model vere conducted by the 
USDA Soil Conservation Service and Resource Consultants, Inc. , in a cooper- 
ative agreement vith NASA. The first part of the report discusses the 
results of application of the model on the Rio Grande vatershed above Del 
Norte in southwestern Colorado. This vatershed encompasses the South Fork 
of the Rio Grande vhich vas studied previously by Shafer et al. (1981a). 

The second part of this report deals vith initial attempts at development 
of a real time predictive mode version of the Martinec-Rango model. Tue 
purpose of this development is to utilize this model to predict mean daily 
flovs due to snovmelt runoff from mountainous vatersheds. Primarily, this 
portion of the study focused on development of methods to predict the 
snov-cover depletion curves for operation of the predictive mode model. 
Recommendations for the implementation of this predictive mode model are 
made. 

The computer listing for the Martinec-Rango model as applied in this report 
is shovn in Appendix A. 
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LARGE BASIN STUDY 

RIO GRANT E RIVER BASIN ABOVE DEL NORTE 

The Rio Grande has its headwaters in the San Juan Mountains of southwest 
Colorado (Fig. 1). This watershed was studied above the stream-gaging 
station Just a few miles west of the town of Del Norte, Colorado. Within 
the boundaries of this watershed is the South Fork of the Rio Grande in 
which the snowmelt runoff was previously simulated using the Martlnec-Rango 
model (Shafer et al., 1981a). 

WATERSHED DESCRIPTION 


The watershed of the Rio Grande River above the Del Norte gage has an area 
of approximately 1,320 square miles. Its elevation ranges from 7,980 feet 
at the stream gage to 13,830 feet at the highest point. Fig. 2 shows the 
limits of the watershed above the gaging station and also shows the eleva- 
tion contours and the division of the watershed into three elevation zones 
chosen roughly to correspond to the mixed conifer-aspen band (Zone A) , the 
spruce-fir band (Zone B) and the area above timberline (Zone C). 

Fig. 3 is a graphical hypsometric analysis of the Rio Grande above Del 
Norte basin. Fig. 4 shows the same information for the South Fork of the 
Rio Grande watershed. Mean elevations of each zone are shown on this 
figure. Table 1 presents a numerical synopsis of the hypsometric curve 
presented in Fig. 3. 


Table 1 

Area-Elevation Statistics 
Rio Grande Watershed above Del Norte 



Rio 

Grande 


Zone 

Elevation Range 
(ft.) 

Area 
(mi J 

X of 

watershed 

A 

7,980 - 9,600 

301 

22.8 

B 

9,6C^ - 11,000 

496 

37.6 

C 

11,000 - 13.380 

523 

39.6 

Total 


1,320 

100 


This watershed is typical of the region where a permanent snowpack begins 
accumulating in late October and generally reaches a seasonal maximum near 
the first of April. During April snowpack depletion is normally exper- 
ienced at lower elevations below 9,500 feet while accumulation may continue 
at middle and high elevations up to the first of May. Winter snowpack 
accumulation is a gradual process which results from many snowfall events 
as contrasted with areas where a few major storms contribute the bulk of 
the permanent snowpack. Snowpacks in both watersheds are sub-freezing 
throughout most of the winter. This heat deficit must be made up and the 
snowpack brought to isothermal conditions in the spring before appreciable 
snowmelt runoff can begin. 
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HYPSOMETRIC ANALYSIS 
SOUTH FORK RIO GRANOE 


MEAN 

ELEVATION 
ZONE C 
11522ft (3512 m)- 


MEAN ELEVATION 
ZONE A 
8996 ft (2742 m) 


8,222 


MEAN ELEVATION 
ZONE B 

10312 ft ( 3143 m) 



AREA (mi 2 ) 


Figure 4 - Hypsometric Analysis for South Fork Rio Grande 
at South Fork, Colorado 
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Snowmelt run f product* 75 to 85 percent of th* Annual flow In tha 

watershed. Thla streamflow fro* waiting anovpacka ia tha lifablood of tha 
Rio Cranda valley, without which atrlcultural production would ba virtually 
nonexistent. It la thla dapandanca on tha enowmelt ru off antaring tha Rio 
Grand* vallay that mekee tha accurat* and timely foracaating of ita 
quantity and temporal diatrlbution ao Important, Critical management 
decisiona allocating the anticipated runoff muat ba mad* before and during 
the major runoff period. Thaaa dacialona are covered by atatutaa and 
regulationa spalled out in state law, interstate, and fntamational 
agreement a. 

The location of pertinent climatological stations, and snow courses in and 
near the study watersheds are shown in Fig. 2. Climatological records for 
each of th stations for the Apr il-Sept ember period is presented In 
Appendix B. 

Streamflow records for th* Dei Norte and South Fork gages adjusted for 
upstream reservoirs ' *prii-September for 1973-80) are given in Appendix C. 
These records show that during this 8-year period wide extremes in runoff 
were observed. A frequency analysis of streamflo-i reveals the drought 
conditions which prevailed during 1977 had a recurrence interval of nearly 
100 years. On the other hand, 1979 was an extremely heavy snowpack year 
which produced runoff that would only be expected to occur once in every 20 
years. Other years fell somewhere between these two extremes. It is 
fortunate for the purposes of this investigation that such diverse hydro- 
logic conditions were experienced in ao short a span of time. Due to thia 
diversity it was possible to exercise the model thoroughly and test its 
suitability in simulating widely disparate snowmelt runoff occurrences. 

Thia watershed contains three reservoirs of hydrologic importance— Rio 
Grande, Santa Maria, and Continental. The effect of the*.* reservoirs on 
the observed flow at th* Del Norte gag# was computed according to analysis 
of reservoir releases and travel times. A short computer program used to 
adjust th# Del Norte observed flow to an adjusted "free natural flow" Is 
also given in Appendix C. 

S TREAMFLOW LAG CHAR ACTERI STICS 

An analysis of the daily hydrographa for thia watershed was made in an 
effort to evaluate the time it took meltwater to reach the South Fork and 
De’ Norte stream gages. Fig. 5 shows a segment of the dally hydrograph 
record for the South Fork Rio Grande near the time of the seasonal peak in 
early June. Fig. 6 shows a similar segment at the end of the snowmelt 
season in early July. During high flow periods approximately /0 percent of 
the flow of day n is a result of melt on the previous day with only about 
30 percent coming from day n's generated melt. These values were obtained 
by planimetering the areas under the daily hydrograph. In July the travel 
time for meltwater to reach the gage is longer due to the recession of the 
snowpack line to a higher elevation. At the lower stage about 80 percent 
of day n's flow ia a result of the previous day’s melt and only 20 percent 
is a result of the current day's melt. At high flows on the Rio Grande 
near Del Norte approximately 65 percent of day n's streamflow was produced 
by day n-l's melt, while at lower flows late in the runoff season, this 
contribution increased to over 90 percent. 
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SOUTH FORK RIO GRANDE AT SOUTH FORK 



Figure 5 - Daily Hydrograph of South Fork Rio Grande at 
South Fork for June 2, 3, 1978, near Peak Period of 
Snowmelt Runoff 


SOUTH FORK RIO GRANDE AT SOUTH FORK 



TIME 


Figure 6 - Daily Hydrograph of South Fork Rio Crande at 
South Fork for July 5, 6, 1978, near End- of Active 
Snowmelt Runoff Recession 
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This gradual shift in the tine distribution of the dally streamflow hydro- 
graph from June to July is a function of the size and shape of the water- 
shed areal distribution of snow cover, and the stage of the stream. In the 
interests of simplicity, it was decided that each day's meltwater produc- 
tion in the South Fork Rio Grande would be split into two components and 
summed to match the time frame of published dally streamflow records, i.e., 
midnight to midnight. The procedure Involved summing 70 percent of melt 
generated on day n-1 and 30 percent of the melt on day n to arrive at the 
total meltwater contribution of day n. Although this technique results in 
the introduction of a small error at low flows it was felt to be insigni- 
ficant enough to ignore. 

On the Rio Grande at Del Norte station, an alternate approach was devised 
wherein a relationship was developed between the preceding day's snowmelt 
contribution to scisamflow (expressed as a percent) and the total basin 
snow cover (also expressed as a percent). This technique takes into 
consideration the size of the basin and the increase in travel time for 
meltwater to reach the gage as the snow line recedes. This relationship is 
graphically illustrated in Fig. 7. This procedure permits a continuous 
relationship rather than the "step-function" approach previously described 
for the South Fork Rio Grande. 

SNOW COVER DEPLETION CURVES 


Basin snow-covered area is a critical variable in the Martinec model. 

Daily values of the proportion of each zone covered by snow are required to 
drive the model. To obtain this information Landsat imagery was analyzed 
by photointerpretive techniques, and a map was constructed delineating the 
location of snow-covered areas in each zone. These areas were then 
manually planimetered and a value obtained for snow cover in each zone. 

All available imagery for each year was mapped and a series of snow-cover 
depletion curves were developed similar to those shown in Fig. 8 for the 
Rio Grande for 1979. Daily values of zonal snow cover are extracted from 
the depletion curve and used as input to the model. Appendix D contains a 
complete summary of snow-cover mapping data and snow-cover depletion curves 
for the South Fork Rio Grande and Rio Grande above Del Norte for the 
1973-80 period. 

A comparison of the snow-cover depletion curves for a given watershed over 
a period of years shows them to be basically similar in shape but shifted 
in time. This time shift is primarily a function of the volume of snow 
stored in the watershed. The relative displacement of the curves from year 
to year are useful in providing a rough measure of the amount of snow on 
the watershed and also the volume of water to be produced when it melts. 

An accurate portrayal of this temporal distribution of the snow-cover 
recession is the most important single factor contributing to the success 
of the model in explaining the variability of snowmelt discharge rates. 

TEMPERATURE LAPSE RATES 


Climatological data for this study was obtained from published records for 
the Del Norte and Wolf Creek IE stations. Since Del Norte (elevation 7,884 
ft.) is located near the lower end of the Rio Grande watershed and the Wolf 
Creek IE station is located on Wolf Creek Pass in the Rio Grande watershed 
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Figure 7 

TOTAL BASIN SNOW COVERED AREA VS. 
PRECEDING DAY'S SNOWMELT CONTRIBUTION 
TO STREAM FLOW 

RIO GRANDE BASIN , 1975 -1979 DATA 
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at an elevation of 10,642 ft., actual daily lapse rates calculated between 
these two stations were utilised to develop degree days for the hypsometric 
mean elevations of each zone. 


To account for variations in station observation time, adjustments were 
made in daily data. Wolf Creek IE observation time is 8 a.m.; therefore, 
daily maximum temperature and precipitation were obtained from the follow- 
ing day's published data. For Del Norte the observation time is 6 p.m.; 
therefore no adjustment was made. 

Temperatures for each zone on the Rio Grande watershed were calculated from 
observed climatological records from Del Norte and Wolf Creek IE according 
to the following equation: 


T « T ♦ ^ DN (T -T) 

ZN DN — U WC DN ; 


where: 


Sic Sn 



temperature for a given zone (°F) 
temperature at Del Norte (°F) 
temperature at Wolf Creek IE (°F) 

hypsometric mean elevation for zone A, B, or C (ft.) 
elevation of Dei Norte station (7,884 ft.) 
elevation of Wolf Creek IE station (10,642 ft.) 


Consideration was given to incorporating the Hermit 7 ESE temperature data 
since it is representative of a rather significant area of the Rio Grande 
watershed. Anlyses showed, however, that this station is greatly influ- 
enced by temperature inversions. This means that it is very difficult to 
use this station for extrapolating temperatures either up or down in 
elevation. 


PRECIPITATION 


Precipitation input for the model was based on data from Wolf Creek IE and 
Del Norte. Precipitation in each zone was determined by two different 
methods. To account for orographic type storms which occur in the winter 
and spring, the precipitation data were calculated by the same elevation 
ratios used for determining lapse rates for temperatures in each zone. 

This method was utilized for precipitation occurring in April and May of 
each year. The following equation summarizes this method: 


+ e zn - e dn 

ZN DN = — vr WC DN ; 


^fC " 


"DN 


where: 

P ZN - precipitation for a given zone (in.) 

P~ N * precipitation at Del Norte (in.) 

P W C * precipitation at Wolf Creek IE (in.) 

E Z n* Ejjjj* E^£ are as defined in the equation just above 
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The second method used to develop precipitation for each zone accounted for 
convection type storms which occur during the late spring and summer. For 
this case the precipitation occurring at Del Norte was utilized directly 
for Zones A and B, and Uolf Creek IE data were used for Zone C. This 
method was used for the months of June through September. In other basins 
alternative precipitation transfer functions may have to be developed to 
generate realistic precipitation amounts in the various zones. 

Appendix A gives the actual published data for Wolf Creek IE and Del Norte 
used in the study. Appendix E contains the actual degree-day and precipi- 
tation data calculated for each zone for the South Fork Rio Grande and Rio 
Grande during 1973-80 snowmelt periods. 

MELT RATES 


Daily snowmelt depths are calculated by the model using the degree-day 
method. Although several approaches in applying this method have been 
advanced by various authors including Gartska (1958), Martinec (1960) and 
Linsley and Franzini (1979), they all relate the amount of snowmelt in a 
basin to a degree-day index factor. The degree-day factor for snowmelt 
computations is normally the positive departure of the mean daily tempera- 
ture above a base temperature of 32 F. The factor is thus an empirical 
measure of the amount of energy available to melt snow based on air temper- 
ature alone. However, other variables including solar radiation, wind and 
humidity also influence melt rates. For this reason a wide variability in 
melt rates based on degree-days can be expected. Melt rates ranging from 
0.05 to 0.20 inches /degree-day (°F) hava been observed at the Upper San 
Juan SN0TEL site (Shafer et al., 1981a). 

In the most common application of the degree-day method, daily mean temper- 
ature is approximated as the average of the observed daily minimum and 
daily maximum temperature for a 24-hour period. This method can lead to 
instances where melt is observed when the number of degree-days is zero. 
Such occurrences are caused by the equal weighting of the dai'y minimum and 
maximum temperatures. 

It is not possible to completely eliminate problems of the foregoing nature 
when attempting to approximate mean daily temperatures from available 
climatological records. Even with its obvious shortcomings, the degree-day 
method does provide a useful tool for modeling. Since it is only an index 
and not an exact relationship, it is not critical in the model. It is 
Important, however, that whatever degree-day method is once chosen, it be 
consistently applied thereafter. 

For the purposes of the current investigation a degree-day was calculated 
on the following basis: 

T + T . - 32°F 

max min 



where: 

T^ * degree day 

T° * maximum daily temperature (°F) 

XT * minimum daily temperature (°F) 
rain 
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This method of calculating degree-days differs from the original model's 
algorithm which uses an effective minimum of 32°F resulting in higher 
degree-days early in the snowmelt season. Estimates for the degree-day 
melt factor (a) were made based on plots similar to Fig. 9 as well as 
snow-course data in each of the various elevatlonal zones of the water- 
sheds. A wide range of melt rates were calulated for different periods in 
different years. A watershed average of about 0.08 inches /degree-day were 
computed for both the South Fork Rio Grande and Rio Grande above Del Norte 
for all years in the simulation. This compares favorably with average 
values found in the literature (U.S. Army Corps of Englneers t 1960). 

A summary of the melt rate factors by zone for each watershed is contained 
in Appendix F with each year's simulation results. 

Melt rates for Zone A in both watersheds generally ranged between .05 
inches/degree-day in April to .07 inches/degree-day for the May-September 
period. Zone B melt rates generally ranged from .03 inches/degree-day in 
April to a maximum of 0.12 lnches/degree-day in late May and early June. 
Zone C melt rates ranged from .02 inches/degree-day in early April to a 
maximum of 0.16 inches/degree-day in early June. 

RUNOFF COEFFICIENT 


A proper evaluation of the zonal runoff coefficient* c, is necessary to 
achieve a reasonable degree of success with the model. An explicit evalua- 
tion of c is not easily performed. Rather, an estimate based on prior 
experience in other watersheds and consistent with hydrologic conditions 
prevailing in the study watershed is normally made. For both the South 
Fork Rio Grande and Rio Grande above Del Norte, estimates of c were 
adjusted between zones and seasonably varied for each year to achieve as 
close a fit as possible between the observed and simulated hydrographs for 
each year. Appendix F contains a listing of runoff coefficients for each 
yearly simulation run for each watershed. 

HYDROGRAPH RECESSION ANALYSIS 


To correctly apportion the observed daily streamflow discharge between the 
current day's snowmelt runoff contribution and the effect of previous day's 
melt requires an evaluation of the recession coefficient, k. Martinec and 
Rango (1979) have shown that this coefficient is related to discharge by 
the following exponential function: 


k 

where: 


» x • Q ^ 

k * recession coefficient 

x and y are parameters unique to a given basin 
Q ■ current discharge in ft /sec 


A direct evaluation of this relationship is possible for the watershed. To 
obtain the values of x and y, values of Q were plotted against Q for 
recession periods over a number of years and a lower envelope line and an 
average line drawn through the points. The mean line is a visual best fit 
approximation. The lower envelope line does not necessarily Include all 
points. It is placed to include approximately 90-95 percent of all points. 
Fig. 10 shows the recession plots for the South Fork Rio Grande and Fig. 11 


15 




INSCHAMC 0-., (eft) 


ORIGINAL PAGE tt| 
OF POOR QUALITY 



Fifur* 10 - Snovaelt Runoff Recession Plot for South Fork Rio Grand* 
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Fig. 11. Snowmelt runoff recession plot for Rio Grande River 
at Del Norte. 
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for Rio Grande near Del Norte. A log-log plot of k versus Q baaed on the 
previous plots was then constructed and the watershed constants k and y 
evaluated for both the envelope line and the mean line. Figures 12 and 13 
show the derived recession relationship as a function of discharge for 
South Fork Rio Grande and Rio Grande near Del Norte, respectively. 

RESULTS OF THE MODEL RUNS 


The model was operated on the Rio Grande in two modes: it was run with and 

without the routing corrections which dictate the amount of the preceding 
day's meltwater contributing to the current day's flow as a function of the 
total snow cover on the entire basin. Comparisons between the computed 
flows and the actual flows with and without the continuous routing function 
are shown for individual years from 1973-79. Those without the correction, 
that is relying on the former step function technique are shown in Figures 
14 to 20, while those with thi corrections are shown in Figures 21 to 27. 
Figures 28 to 34 show Actual md simulated flows on the South Fork using 
the step function only. 


Numerical analyses of these model runs have been made using the Nash- 
Sutcliffe correlation technique. The Nash-Sutcliffe (1970) R is somewhat 
analogous to the coefficient of determination and is defined as follows: 


n E ( V q) - 
i - 1 


n E 
i 


(q,- q< x ) 


I ? <V fl > 2 
n i - 1 

where ^ 

R is a measure of model efficiency 
q^ ■ observed discharge 

^i ■ simulated discharge 
q - mean of observed discharge 
n - number of discharge values 

This technique, using in time series correlations, tends to give greater 
emphasis to the higher values. With this technique it is possible to 
obtain a coefficient of determination less than zero. These results are 
presented in Table 2 along with an analysis of the percentage seasonal 

difference ( — c tua l " - Computed . * ^ positive value Indicates an over- 
prediction. actual 

Analysis of Figures 14-27 and results shown in Table 2 reveal nearly 
identical performance is achieved by the model with either the step funC' 
tion or continuous function routing. As a practical matter, it would 
therefore be preferable from an operational perspective to use the con- 
tinuous function to minimize parameter entry. 

Table 2 also includes a summary of the model's performance for the South 
Fork Rio Grande using the step function routing only. 
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Figure 12 - Snowmelt Runoff Recession Coefficient for South Fork Rio Grano* 
as a Function of Discharge 
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Fig. 13, Snowmelt coefficient for the Rio Grande River at Del Norte 
as a function of discharge. 
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RIO GRANDE RIVER NEAR OELN0RTE, 1973. 

WOLFCR IE AND DELNORTE DATA WITHOUT ROUTING CORRECTION. 
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Fig. 15 

RIO GRANDE RIVER NEAR DELN0RTE. 1974. 

H0LFCR IE AND DELNORTE OATA H1THOUT ROUTING CORRECTION. 
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HOLFCR IE AND OELNORTE DATA WITHOUT ROUTING CORRECTION 
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Fig. 17 

RIO GRANDE RIVER NEAR DELNORTE. 1976 

WOLECR IE AND DELNORTl DATA WITHOUT ROUTING CORRECTION 
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Fig. 18 

RIO GRANDE RIVER NEAR DELNORTE , 1977. 

WOLF CREEK IE AND DELNORTE DATA WITHOUT ROUTING CORRECTION 
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Fig. 19 

RIO GRRNDE RIVER NERR DELNORTE, 1978. 

WOLF CREEK IE AND OELNORTE DRTfl WITHOUT ROUTING CORRECTION. 
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RIO GRANDE RIVER NEAR DELNORTE, 1979. 

WOLF CREEK IE AND DELNORTE DATA WITHOUT ROUTING CORRECTION. 
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Fig. 21 

RIO GRANDE RIVER NEAR 0ELN0RT' . 1973. 

H0LFCR IE AND DELNORTE DATA HITH ROUTING CORRECTION. 
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Fig. 22 

RIO GRRNOE RIVER NEAR DELNORTE, 1974. 

WOLFCR IE AND DELNORTE DATA WITH ROUTING CORRECTION. 
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APRIL HRT JURE JULT RUGUST SEPTEMBER 

Fig. 23 

RIO GRANDE RIVER NEAR DEIN0RTE, 1975. 

HOLFCR IE AND DELNORTE DATA WITH ROUTING CORRECTION. 
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F±g. 25 

RIO GRANDE RIVER NEAR DELN0RTE, 1977. 

WOLFCR IE AND OELNQRTE DATA WITH ROUTING CORRECTION. 
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RIO GRANDE RIVER NEAR DELN0RTE. 1978. 

HOLFCR IE AND DELNORTE DATA WITH ROUTING CORRECTION. 
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Fig. 27 

RIO GRflNOE RIVER NEAR OELNORTE , 1979. 

W0LFCR IE AND OELNORTE DATA WITH ROUTING CORRECTION. 
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SOUTH FORK OF THE RIO GRRNDE RIVER AT SOUTH FORK 
MODEL RUN FOR 1974 RUN MADE OCTOBER 24. 1980 
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SOUTH FORK OF THE RIO UKHNUt rwvcn m suum . u.. 
MODEL RUN FOR 1976 RUN MADE OCTOBER 24, 1980 
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APRIL MAY JUNE JULY AUGUST SEPTEMBER 

FIGURE 32 

SOUTH FORK OF THE RIO GRANDE RIVER AT SOUTH FORK 
MODEL RUN FOR 1977 RUN MADE OCTOBER 24, 1980 
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SOUTH FORK OF THE RIO GRANDE RIVER AT SOUTH FORK 
MODEL RUN FOR 1979 RUN MODE OCTOBER 24, 1980 



Table 2 

Summary of Simulation Results for 
Rio Grande and South Fork 
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PREDICTIVE MODE DEVELOPMENT USING THE 
MARTINEC-RANGO SNOWMELT RUNOFF SIMULATION MODEL 


INTRODUCTION 


Previous work with the Martlnec-Rango model has been restricted to simulat- 
ing streamflow using actual records of climatological and watershed data. 

This phase of the study focused on the ability to use the model in a real- 
time predictive mode. The greatest benefit of any snowmelt simulation 
model is its real-time predictive mode, hence its eventual application in a 
forecast rather than simulation mode. In the predictive mode this model is • 

potentially capable of making short-term projections of the snowmelt hydro- 
graph on the order of several days. This capability is of paramount impor- 
tance for such applications as irrigation water management, reservoir 
management, flood forecasting, and hydropower generation. 

In order to run in the predictive mode, estimations must be made of the 
model parameters and predictions of temperature and precipitation within 
the basin must be obtained from weather forecasts. One key parameter that 
must be estimated for the model is the percent of basin snow cover within 
each elevation zone. Previous studies have used snow-cover data from 
satellite imagery (Shafer et al. 1981b) (Rango & Martinec, 1979), but 
real-time availability of this imagery limits its use when the model is 
applied in the predictive mode. Therefore, other methods must be employed 
to estimate the snow cover. Another parameter necessary for the model is 
the melt rate factor, but, as shown in Shafer et al. (1981a), these can be 
obtained from snow pillow data. The remaining parameters necessary for 
operation of the model in the predictive mode are the runoff coefficients 
for each zone. Model performance is very sensitive to the runoff 
coefficients. Unfortunately, an objective technique to reliably estimate 
this parameter has not been developed in this study. However, it is 
possible to develop a "feel" for realistic values from prior simulations 
and knowledge of the basin's hydrologic response. Additional effort should 
be devoted to development of objective functions to estimate the rt>' off 
coefficients. In the interim, an iteration process over a limited range of 
c values will yield sufficiently accurate estimates to employ in the model. 

Given that all parameters could be estimated and weather forecasts were 
available, it then requires a method of predicting the snow-cover depletion 
curves based on data available in a real-time mode. The objective of this 
effort is to make short period forecasts of the snow cover on each zone 
within the basin. Martinec & Rango (1982) presented methodologies for 
developing snow-cover depletion curves based on cumulative degree days. 

Rango presented families of curves depicting the relationship of snow cover 
to cumulative degree-days in low, average, and high snowpack years. Selec- 
tion of which type curve to use for operational purposes was dependent to a 
large extent on the hydrologist's judgment after weighing several factors 
including the mean basin snowpack water equivalent, the amount of precipi- 
tation occurring on the basin, and the amount of melting that has already 
occurred. Alternative means were sought to more objectively determine the 
snow-cover depletion in widely varying snowpack environments where 
satellite data were lacking. This study investigated several methods of 
predicting the snow-cover depletion curves based on real-time information. 
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ESTIMATION TECHNIQUES FOR 
SNOW-COVER DEPLETION CURVES 


Three new techniques were developed and evaluated for estimating snow-cover 
depletion curves. In order to evaluate the sensitivity of the model to 
these techniques, the information developed In Shafer et al. (1981a) for 
the South Fork of the Rio Grande was used for 1973 through 1979. In 
addition, data present in the previous section of this report for the Rio 
Grande above Del Norte for 1973 through 1979 were also used. 

Multiple Regression Analysis 

The first method examined for predicting the snow-cover depletion curves 
utilized a multiple-regression analysis of April 1 snow-course measurements 
and cumulative degree days in each zone. The degree days were accumulated 
from April 1 for each year, 1973 through 1979, for each zone. Using the 
snow-course data within the basin, estimates of average snow-water equiva- 
lents for each zone on April 1 were calculated and used as the second 
parameter in the analysis. Table 3 shows the snow course used for each 
zone to determine April 1 average water equivalent. Elevational averages 
were computed for the mean water equivalent in each zone. Tne square and 
cube of the accumulated degree days and water equivalents were used in the 
regression analysis. Table 4 presents the results of the regression 
analysis for each zone. Similarly, Table 5 shows the corresponding results 
for the Rio Grande above Del Norte basin. 

Using these predictive regression equations, the model was run for the 1973 
through 1979 period and compared to the results of the actual calibrated 
runs. These results are summarized £n Tables 6 and 7. As can be seen from 
Tables 6 and 7, the Nash Sutcliffe R s decreased for all yearn with the 
poorest fits occurring in 1975 and 1977. In addition, the total predicted 
volume of runoff ranged from a -17 percent to a +18 percent wtich is 
considerably worse than the existing calibrated runs. Appendix G contains 
the plotted hydrographs for these estimates. For both watersheds Figures 
35 through 43 present the actual snow-cover depletion curves for each zone 
for the years 1977, 1978, and 1979 on the South Fork drainage. Figures 44 
through 52 show similar data for the Rio Grande above Del Nort** , The 
estimated snow-cover depletion curves using this technique and the other 
two techniques presented later in the report are also shown. 

Table 3 

Snow-Course Data Used for April 1 Water 
Equivalent in Linear Estimation Technique 


Zone 

Snow 

Course Used^ 

A 

Big Meadow 

Pass Creek 

B 

Upper San Juan 

Wolf Creek Pass 

C 

Wolf Creek Summit 

Highway 


1 / 


Zone water equivalent was estimated by weighted elevation average of 
snow-course data. 
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Table 4 

Multiple Regression Equations Developed 
for South Fork Watershed 

SC - A ♦ B(ACDD) + C(ACDD) 2 ♦ D(ACDD) 3 ♦ E(WE) F(WE) 2 + G(WE) 3 

where ACDD » accumulated degree days from April 1 

WE * water equivalent for April 1 obtained from snow courses lr 
Table 3 

SC * Zone snow-cover percentage 


Constant 

Zone A 

Zone B 

Zone C 

A 

B 

C 

D 

E 

F 

G 

0.2567 
-0.0025190 
-1.40782 x 10"° ( 
1.599848 x 10‘ 
-0.^2629 
0.009516 
-0.00025244 

-3.3705 

-0.0037425 , 

„ 4.9397 x 10"® 

* -2.63288 x 10"* 

0.35312 
-0.0095234 
0.000085909 

0.47112 

-0.00371108 , 

4.6164 x 10 q 
-2.04708 x 10 
0.021492 
-0.0003243 , 
2.468 x 10"® 

R 2 

0.89 

0.88 

0.92 


Table 5 

Multiple Regression Equations Developed 



for Rio 

Grande Watershed 



s.c. - + c 2 (acdd) + 

c 3 (acdd) 2 + c 4 (acdd) 3 + 

c 5 (we) + c 6 (we) 2 


+ C ? (WE) 3 



Constant 

Zone 1 

Zone 2 

Zone 3 

C C 1 

c 2 

C 3 

C 4 

£ 

C* 

/ 

0.24605 
-0.003332 , 

2.22 x 10‘° 

5.571731 x 10 
-0.07765 
0.014701 
-0.00037289 

0.71580 

-0.0047133 . 

7.633 x 10"® 

-10 -9 

-5.18099 x 10 * 

-0.03719 

0.0019075 

-0.00001951 

0.22339 
-0.00477972 , 

7.4739 x 10 

-3.86143 x 10" 9 
0.046839 
-0.00098406 , 
7.163 x 10 

R 2 

0.89 

0.88 

0.94 
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Table 6 

Comparison of Estimation Techniques for Snow-Cover Depletion Curves 
for South Fork of the Rio Grande 
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Shafer et al. (1981a) 


Linear Estimation 


The second technique analyzed for estimation of snow-cover depletion curves 
was based on estimating three points of the snow-cover depletion curve. 
These points were the break point (the point at which the snow cover in the 
zone becomes less than 100 percent), the 50 percent snow-cov*r point, and 
the 0 percent snow-cover point. The linear estimation equations were 
developed using 1973-81 data. 

Predictions of these three points on the curve were made utilizing cumula- 
tive degree days and the same mean zonal April 1 water equivalents of 
snow-course data. Table 3 summarizes the snow courses used for the April 1 
water equivalent estimate. These snow-course values were adjusted by 
altitude to the mean altitude in each zone. Degree days were accumulated 
from the break-point date for estimation of the 50 and 0 percent snow-cover 
dates. The break-point date was estimated using regression analysis of the 
April 1 water equivalents. Summary of the regression equations in this 
technique are shown in Table 8. The accompanying data for the Rio Grande 
above Del Norte is shown in Table 9. 

Intermediate values on the snow-cover curve were then estimated utilizing 
linear relationships between snow-cover depletion and cumulative degree 
days from April 1. The predicted snow-cover curves using this technique 
are shown in Figures 35 through 52 for the years 1977 through 1979. As 
evidenced from these figures the estimation of snow cover is generally 
poorer than utilizing the multiple regression technique. Results of the 
model runs are presented in Tables 6 and 7 and again show similar variation 
with the multiple regression technique. 
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Fig. 35. Actual and Estimated Snow-Cover Depletion Curve — South Fork Rio Grande 
1977 - Zone A. 
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MARCH APRIL MAY JUNE JULY AUGUST 

Fig. 36. Actual and Estimated Snow-Cover Depletion Curve — South Fork Rio Grande, 
1977 - Zone B. 
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Fig. 37. Actual and Estimated Snow-Cover Depletion Curve — South Fork Rio Grande 
1977 - Zone C. 
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Fig. 41. Actual and Estimated Snow-Cover Depletion Curve — South Fork Rio Grande, 
1979 - Zone A. 
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Fig. 42. Actual and Estimated Snow-Cover Depletion Curve — South Fork Rio Grande 
1979 - Zone B. 
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MARCH APRIL MAY JUNE JULY AUGUST 

Fig. 43. Actual and Estimated Snow-Cover Depletion Curve-South Fork Mo Grande. 
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Fi8 ‘ ?Q?7 al and Estimated Snow-Cover Depletion Curve— Rio Grande above Del Norte 

I7/ / - Zone A, 
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MARCH APRIL MAY JUNE JULY AUGUST 

Fig. 46. Actual and F.stimated Snow-Cover Depletion Curve — Rio Grande above Cel Norte, 
1977 - Zone C. 
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Fig. 50. Actual and Estimated Snow-Cover Depletion Curve— -Rio Grande above Del Norte 
1979 - Zone A. 
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Fig. 52. Actual and Estimated Snow-Cover Depletion Curve — Rio Grande above Del Norte 
1979 - Zone C. 
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Table 8 

Summary of Regression Equations Used 
In Linear Estimation Technique for South Fork 

BP « A + B(WE) 

- C + D(WE) 

* E + F(WE) 

where BP ■ Break point day from April 1 (first day of less than 
100% snow cover) 

*7E - Water equivalent on April 1 

= Accumulative degree days from break point to 50% 
snow cover 

(ACDD ) * Accumulative degree days from 50% snow cover to 0% 
° snow cover 


(ACDD_ n ) 

(ACDtfp 




Break point 

(acdd 5Q ) 

(ACDD q ) 

Zone 

A 

B 

2 

r 

C 

D 

2 

r 

E 

F 

2 

r 

A 

-25.83 

1.911 

0.524 

-2.59 

10.811 

0.816 

-2.69 

24.45 

0.511 

B 

-20.66 

1.077 

0.720 

-129.3 

7.264 

0.540 

-9.046 

12.55 

0.681 

C 

-15.65 

1.157 

0.806 

-21.94 

3.765 

0.689 

-85.53 

13.43 

0.806 
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Table 9 

Summary of Regression Equations Used 
in Linear Estimation Technique for Rio Grande 

BP » A ♦ B(WE) 

- C ♦ D(WE) 

- E + F(WE) 

where BP * Break point day from April 1 (first day of less than 
100Z snow cover) 

WE ■ Water equivalent on April 1 

(ACDD^q) ■ Accumulative degree days from break point to 50Z 
snow cover 

(ACDD ) * Accumulative degree days from 50Z snow cover to OZ 


(ACDD sn ) 

(ACDlJJj) 


Break point 

(acdd 50 ) 

(ACDD ) 
o 

Zone A 

B 

2 

r 

C 

D 

m 

E 

F 

2 

r 

A -36.31 

1.781 

0.516 

-21.99 

8.521 

0.713 

180.92 

7.50 

0.236 

B -43.10 

1.399 

0.704 

-35.87 

3.768 

0.641 

118.19 

4.36 

0.309 

C -26.89 

1.037 

0.734 

1.17 

2.310 

0.572 

86.38 

8.15 

0.642 
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PARABOLIC AND EXPONENTIAL ESTIMATION 


The third technique analyzed is similar to the linear estimation technique 
presented previously. The following modifications were made from the 
linear estimation technique: (1) substitution of April-September stream- 

flow volume for water equivalent as a measure of water stored in the form 
of snow; (2) modification of the accumulation periods for degree days; and 
(3) development of a parabolic function for 100 percent to SO percent snow 
cover and an exponential function from 50 percent to 0 percent snow cover 
(hence the notation parabolic and exponential estimation). 

The first modification involved developing a relationship between April- 
September streamflcw volumes and the date of the 50 percent snow cover and 
the date of a 0 percent snow cover. In an operational forecast mode the 50 
percent snow-cover date would be generated from April-September streamflow 
forecasts on April 1. The 0 percent snow-cover date would be determined 
operationally using May 1 seasonal volume forecasts. 

The final modification of the linear estimation technique was to utilize a 
parabolic fit of the accumulated degree days to snow-cover curve from the 
100 percent cover to 50 percent cover using actual data for this interval. 
The snow-cover depletion curve for less than 50 percent cover was predicted 
using an exponential fit relating accumulative degree days to snow-cover 
depletion. This fit was obtained using the observed data in the interval 
of 50 percent snow cover to 0 percent snow cover. Because of the separate 
fitting procedures, a discontinuity exists at the 50 percent snow-cover 
point. These fitting procedures used 1973-81 data. Table 10 summarizes 
the equations used for the prediction of snow-cover depletion curves for 
South Fork, and Table 11 shows similar information for the Rio Grande. 

The results of these estimates are shown in Figures 35 through 52. The 
results of model runs are shown in Tables 6 and 7. 

Generally the results utilizing this method are as good as the results 
utilizing the multiple regression technique. There ar«. .wo advantages of 
this method over the multiple regression technique., The first is the 
ability to change the lower portion of the curve (ACDD ) as a result of 
May 1 estimates of watershed yield, and the second is the ability to adjust 
the estimated ACDD^q and ACDD q with actual imagery during the course of the 
melt season. This adjustment is a result of the normalized equations in 
which snow cover is determined by a ratio of accumulated degree days to the 
accumulated degree days at 50 percent or 0 percent snow cover. If the snow 
cover is known, and the accumulated degree days are known, then the ACDD^q 
and ACDD can be recomputed through the relationships of snow cover and 
accumulated degree days. 
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Table 10 

Summary of Equations Used In 

Parabolic and Exponential Estimation Technique for South Fork 
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OPERATION OF MODEL IN THE 
PREDICTIVE MODE FOR 1980 


To adequately assess the techniques for predicting sncw-cover depletion 
curves* the predictive mode model was run for an actual case study on two 
watersheds. The watersheds selected were the South Fork of the Rio Grande 
at South Fork and the Rio Grande River near Del Norte. These watersheds 
were simulated in the predictive mode for the 1980 snowmelt runoff season. 

The simulations used two techniques, for estimation of the snow-cover 
depletion curves. The first technique utilized the subjective approach 
using a family of snow-cover depletion type curves for each zone (Rango ® 
Martlnec 1982). The second technique was to utilize the parabolic and 
exponential estimation technique to predict the snow-cover depletion curves 
as presented previously. In both cases simplifying assumptions were made 
in order to accurately assess the predictions made by the two techniques. 
The first simplification was to use the actual temperature and precipita- 
tion data for each zone. The second simplification was to assume that 
satellite imagery would be immediately available to update the predictions 
of snow cover. 

The predictions of snow cover were updated at the dates of actual satellite 
imagery. For the case of the subjective technique, the type curve merely 
is adjusted for predictions of further snow-cover depletion data until the 
date of the next image. For the case of the estimation technique, the 
snow-cover estimates obtained from the parabolic and exponential regression 
equations were updated with data from actual imagery to revise the estimate 
of the accumulative degree days to 50 percent or 0 percent, respectively. 

In order to minimize the effects of parameter estimation for this study, 
the 1980 actual snow-cover data were utilized for both watersheds in order 
to determine snowmelt rate factors and runoff coefficients (Figures 53 and 
54). Results of these so-called "actual" simulations are shown in Table 12 
and are plotted in Figures 55 and 56. 


TYPE CURVE TECHNIQUE 

This procedure involve selecting a "type" curve from among a family of 
modified depletion curves relating snow-covered area to accumulated degree 
days starting April 1. Fig. 57 is an example of such a family generated 
for zone A of the South Fork of Rio Grande. In an operational mode one 
would select the type curve based upon knowledge of the existing snowpack 
conditions and the projected seasonal volume runoff. For example, if a 
much above average snowpack existed on April 1 and seasonal streamflow 
forecasts indicate an abnormally high runoff, one would select a type curve 
which was developed from a similar magnitude year in the past. As satel- 
lite imagery becomes available during the course of the snowmelt period, 
adjustments can be made in the shape of the type curve using actual temper- 
ature data. This approach is simplistic in that it presupposes the same 
precipitation amounts in the current year as occurred in the year chosen as 
the appropriate type curve. Thus, any deviation in precipitation pattern 
compared to the year selected will be a source of error. 
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To test the reliability of this method, predicted snow-cover depletion 
curves were constructed for 1980 for the South Fork Rio Grande and Rio 
Grande above Del Norte using actual temperature data and families of 
modified depletion curves based on the 1973-1979 period. As satellite 
Imagery became available, the predicted snow-cover depletion curves were 
adjusted to conform to actual values. The predicted snow-cover depletion 
curves were then used as input data to the model. Figures S3 and 54 show 
the actual and predicted snow-cover depletion curves for the South Fork Rio 
Grande and Rio Grande, respectively. 

Table 12 shows the results of the model runs using this technique. The 
computed hydrographs for South Fork and Rio Grande are shown on Figures 58 
and 59, respectively. 

Table 12 

Results of Predictive Mode Runs for 1980 


South Fork 
N-S R Z 

Run 1/ 

of Rio Grande 
Seasonal 
difference 
(X) 

Rio Grande 
N-S R^ 

U 

above Del Norte 
Seasonal 
difference 
(X) 

"Actual" simulation 

0.9177 

0.27 

0.8978 

0.44 

Type curve technique 

0.8344 

14.69 

0.5656 

23.56 

Parabolic and 
exponential technique 
without satellite 
update 

0.7719 

3.47 

0.8674 

1.97 

Parabolic and 
exponential technique 
with satellite update 

0.2359 

24.87 

0.4872 

21.74 


17 2 

— ' Nash-Sutclif fe R 
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Fig. 53. Actual and Predicted Snow-Cover Depletion Curve, 
South Fork Rio Grande - 1980. 
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Fig. 54. Actual and Predicted Sno'*-Cover Depletion Curve, 
Rio Grande above Del Not te - 1980. 
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Fig. 55. 

SOUTH FORK OF THE RIO GRRNOE. RT SOUTH FORK. 1980. 
ACTUAL SNOW COVER OR T A. 
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SOUTH FORK OF THE RIO GRANDE , AT SOUTH FORK. 1980. 
TYPE CURVE METHOD USED FOR ESTIMATE OF SNOW COVER. 
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Fig. 59. 

RIO GRRNOE RIVER. NEAR DEL NORTE. 1980. 

TYPE CURVE METHOD USED FOR ESTIMATE OF SNOW COVER 


PARABOLIC AND EXPONENTIAL TECHNIQUE 


The parabolic and exponential technique discussed previously was used to 
predict the snow-cover depletion curves for both the South Fork and Rio 
Grande watersheds. The equations presented previously were utilized to 
calculate the snow-cover depletion based on April 1 and May 1 snow survey 
forecasts of watershed yields and using actual degree day information lor 
each zone. 

The parabolic and exponential equations were solved in reverse on dates 
where satellite imagery was available to update snow-cover information. 

The equations were solved for revised estimates of accumulated degree days 
to the 50 percent and 0 percent snow covers. Several assumptions were made 
in order to make the revisions to the curves. The first assumption was to 
adjust the parameters to allow for a 95 percent snow cover when satellite 
imagery showed there was 100 percent snow cover on a zone. If 100 percent 
was used in the equations, a value of inlinity would be obtained for the 
ACDD^q parameter. 

The advantage of this technique over the type curve technique is the 
ability to allow the computer model to generate its own snow-cover data 
base upon input of actual measured data. This is desirable in the 
operation* 1 application of the model provided the algorithms developed 
produce satisfactory estimates of the streamflow hydrograph. 

The results of the runs for this technique, the type curve technique, and 
the actual simulation technique (utilizing actual s..ow-cover data) are 
summarized in Table 12 for South Fork and Rio Grande, respectively. Plots 
of actual and estimated snow-cover depletion curves are shown in Figures 53 
and 54. 

Computer plots of each run summarized in the tables are shown in Figures 60 
through 63. 

Runs made utilizing the parabolic and exponential technique without satel- 
lite imagery revision provided better fits of the actual streamflow hydro- 
graphs for both watersheds than when the satellite imagery was used to 
revise the estimates. Further investigation of this approach will be 
necessary to improve upon the updating procedure. 

Results cf this analysis indicate that for the single test year, the type 
curve technique provided a better fit of the actual seasonal hydrographs. 

The parabolic and exponential method was superior when no satellite imagery 
was used but was less suitable when updated with actual satellite snow- 
cover data. This indicates that the algorithm used for revision of snow- 
cover depletion curves for the parabolic and exponential t^hnique must be 
revised to better update the estimated curves with acti 1 . atellite imagery. 
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S0UTH°F0RK OF THE RIO GRRNDE RIVER, RT SOUTH FORK. 1980. 
PARBOLIC AND EXPONENTIAL ESTIMATE OF SNOW COVER. 
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RIO GRANDE RIVER NEAR DEL NORTE . COLORADO. 1980. 

PAR. AND EXP. WITH SATELITE UPDATE ESTIMATE OF SNOW COVER. 


SUMMARY 


The Martinec-Rango snowmelt runoff model was applied to two watersheds in 
the Rio Grande basin, Colorado. The first watershed was the South Fork Rio 
Grande, a medium size drainage encompassing 216 ml without major reservoirs 
or diversions. The second watershed^was the Rio Grande above Del Norte, a 
large drainage encompassing 1,320 mi and containing three major reservoirs. 
The model was successfully applied to both watersheds when run in a simula- 
tion mode for the period 1973-79. This period included both high and low 
runoff seasons and was a good test interval. 

The model performed better on the South Fork than on the Rio Grande main 
stem as expected for two reasons: (1) the streamflow records on the Rio 

Grande main stem were adjusted for upstream reservoir regulation which 
induced an added source of error, and (2) the climatic data stations used 
for both watersheds were wholly in, or in closest proximity to, the South 
Fork drainage. It is believed that had additional climatic ac^a been 
available for the upper reaches of the Rio Grande watershed, the simulation 
results on both watersheds would have been similar. The average variance 
accounted for by the model in the 1973-79 period was 86.2 percent for the 
Rio Grande watershed and 89.0 percent for the South Fork. These results 
demonstrate the inherent flexibility in the model for application to 
watersheds of widely disparate size. This characteristic is highly desir- 
able because it enhances the model's capacity for universal application. 

A second aspect of the study focused on adapting the model to rui. in a 
forecast mode. Central to this effort was the need co develop a technique 
to forecast the shape of the snow-cover depletion curves between satellite 
passes or in the absence of snow-cover information for an extended period. 
Four separate approaches were investigated. They include simple linear 
estimation, multiple regression, parabolic-exponential, and type curve. 
Preliminary analysis indicated only two of the four approaches considered 
were worthwhile candidates for a trial application in an operational 
context. Both the parabolic-exponential and type curve methods were run on 
the South Fork and Pio Grande watersheds for the 1980 runoff season using 
satellite snow-covcr updates when available. 

The results of the trial were disappointing. Neither method provided an 
acceptable fit for operational forecasts. This result may be partially 
attributable to a large amount of estimated climatic data for one of two 
primary base stations during the 1980 season. However, it is evident from 
the analysis that none of the four techniques in their present form are 
entirely satisfactory for predicting snow-cover depletion operationally. 
Further investigations are recommended before adoption of any standardized 
snow-cover prediction technique. 

Preparatory to using the model in a forecast mode a routine for updating 
with observed streamflow is essential. Although no effort was made in this 
study to include this provision it would be a relatively simple task to 
accomplish and would increase the model's utility. 
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Computer Model Listing 


FORTRAN IV 


V02.1-1 


Tue 25-Ha*-82 16J46:i4 


FACE 001 
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0011 
0012 
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c 

c 

c 

c 

c 
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PROGRAM SNOMEL 


THIS PROGRAM PREDICTS DAILY STREAMFLOW FROM TEMFERATURE 
AND SNOWCOVER DATA FOR A GIVEN WATERSHED AREA 

THE WATERSHED SHOULD DE DIVIDED INTO ELEVATION 
ZONES. THIS VERSION CAN ACCEPT UP TO FIVE ZONES 

THE PROGRAM SIMULATES THE SIX MONTH PERIOD F.>OM 
APRIL 1 THROUGH SEPTEMBER 30. ALL INPUT DATA 
IS REQUIRED FOR THIS PERIOD ONLY. 

THE VARIABLES USED ARE • 

T = TEMFERATURE IN DEGREE F/DAY ABOVE 32 DEG F. 

DEVELOPED FOR EACH ELEVATION ZONE 
P = THE PRECIPITATION IN INCHES/DAY FOR EACH ZONE 
S = THE SNOW COVER AREA IN EACH ZONE IN % OF 
THE TOTAL AREA OF THE ZONE 

A = THE SNOWMELT CGEFICIENT IN INCHES/DAY PER DEGREE F 

C = THE RUNOFF COEFICIENT FOR EACH ZONE 

AREA = THE AREA OF EACH ZONE IN SQ. MILES 

STRFLO = THE COMPUTED DAILY STREAM FLOW 

ACTFLO = THE ACTUAL STREAM FLOW MEASURED 

YEAR = THE YEAR BEING MODELED 

NZ = THE NUMBER OF ELEVATION ZONES 

DATA REQUIREMENTS ARE.* 

THE DAILY TEMPERATURE IN EACH ZONE IN DEGREE DAYS 
THE DAILY PRECIPITATION DATA FOR EACH ZONE 
THE DAILY SNOW COVER DATA FOR EACH ZONE 
THE SNOWMELT AND RUNOFF COEKICIENTS FOR 
15 DAY INCREMENTS BEGINNING APRIL 1 
NUMBER OF ZONES AND THE AREA OF EACH ZONE 
RECESSION COEFICIENTS FOR THE GAGING STATION 


THIS MODEL WAS DEVELOPED BASED ON WORK BY RANGO AND MARTINEC 
DEVELOPED MARCH 1980 BY RESOURCE CONSULTANTS INC. 

VERSION 1.4 - REVISIONS OF OCTOBER 1981 


DIMENSION T<5> 183) >S<5» 183) »F<5» 183) r A<5, 1 83 ) t C < 5> 1 83 ) » AREA ( 5 ) 
DIMENSION RUNOFF <5 1 183 ) t STRFLO < 1 83 ) > ACTFLO < 183 ) , FT; ( 5 . 1 83 ) 
DIMENSION NDAY (6) > NAME < 20 ) > DATE < 3 ) 

DIMENSION M0<12>3) 

DOUBLE PRECISION RECESS » F'RUNOF » RUNOFF > QN iX»Y» SUMRO 
INTEGER*4 MO > YEAR t NAME * DATE » PR 
DATA M0( 1 » 1 ) » M0< 1 » 2 ) » M0< 1 » 3 ) /4H 
M0<2» 1 ) >M0<2»2) »M0<2t3)/4H 
MO ( 3 » 1 ) » MO < 3 * 2 ) f MO <3>3)/4H 
MQ< 4 » 1 > »M0<4*2) »M0(4.3)/4H 
M0( 5 * 1 ) » MO ( 5 » 2 ) t M0< 5 » 3 ) /4H 


DATA 

DATA 

DATA 

DATA 

DATA 


MO ( 6 ? 1 ) »M0<6»2) »M0(6.3)/4H 


JA * 4HNUAR f 4HY 
FE » 4HBRUA » 4HRY 
MA » 4HRCH t 4H 
AP t 4HUIL » 4H 
MA » 4HY t 4H 

JU » / FINE t 4H 


/ 

/ 

/ 

/ 

/ 

/ 


C " 
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Tue 25-M**-82 16.‘46d4 


PAGE 002 


FORTRAN 10 


002 . 1-1 


0014 

0015 

0016 

0017 

0018 
001 ? 
0020 


0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 


DATA MO < 7 * 1 ) » NO (7.2). MO (7.3) /4H JU.4HLY »4H / 

DATA M0(8» 1 ) >M0<8»2) »M0(8»3)/4H AU» 4HGUST . 4M / 

DATA M0< 9* 1 ) »M0( 9. 2) »H0( 9. 3)/4H SE » 4HPTEM » 4HBER / 
DATA M0( 10t 1 > »M0( 10 » 2) >M0< 10»3)/4H 0» 4KCT0P » 4HER / 

DATA MO (11.1). MO (11.2). MO ( 1 1 • 3 ) /4H NO » 4HVEMB » 4HER / 
DATA M0< 12* 1 > f HOC 12 » 2) »M0( 12» 3)/4H D£ » 4MCEMB » 4HER / 
DATA NDAY/30.31 »30»31 .31 *30/ 


C 

C READ IN THE WATERSHED PARAMETERS 

C...... CARD i: 

C...... NAME OF WATERSHED (NAME) 20A4 COLUMNS 1-80 

C CARD 21 

C NUMBER OF ZONES (NZ) 12 COLUMNS 1-2 

C AREA OF EACH ZONE (AREA) 5F10.3 COLS 3-52 

C CARD GROUP 3J 

C SNOWMELT COEF. (A) 12F6.4 COLUMNS 1-72 

C ONE CARD FOR EACH ZONE <2 VALUES PER MONTH) 

C CARD GROUP 4 i 

C RUNOFF COEF. (C) 12F6.4 COLUMNS 1-72 

C ONE CARD FOR EACH ZONE (2 VALUES PER MONTH) 

C CARD GROUP 5J 

C PREC. RUNOFF METHOD (PR) 611 COLUMNS 1-12 

C PR=OJ PREC. RUNOFF IS COMPUTED FROM 

C NON-SNOWCOVERED AREAS ONLY 

C PR=1J PREC. RUNOFF IS COMPUTED FROM 

C THE TOTAL ZONE AREA 

C ONE CARD FOR EACH ZONE (2 VALUES PER MONTH) 

C CARD 6i 

C RECESSION COEF. (X.Y) 2F10.5 COLUMNS 1-20 

C DAILY RUNOFF REACHING GAGE (DR = RA + RB(SNOWCOVER 

C 2F10.5 COLUMNS 21-40? 

C 


CALL ASSI6N( 1 » ' WATSHE . DAT ' . 10* ' RDO ' ) 

REAIKl.lOl) NAME 

101 FORMAT (20A4) 

READ( 1 » 102) NZ » AREA 

102 FORMAT ( I2»5F10.2) 

DO 1 J= l.NZ 

1 READ(1 » 103) A(J»1)»A(J»16)»A(J»31)»A(J»46)»A(J» 62 >»A(J»77).A(J» 92 
1»A(J»107)»A(J»1 23 ) » A < J * 1 38) »A(J»154) »A(J»169) 

DO 5 J=1.NZ 

READdf 103) C(J»l)fC(J»16)»C(J»31)»C(J»46)»C(J»62)»C(J»77)»C(J» 92 
1*C( J»107> »C< J>123) »C< J» 138) vC< J»154> * C< Jf 169) 

5 CONTINUE 


DO 40 J=1 »NZ 

READ(1» 116) PR ( J » 1 ) » PR ( J 1 1 6 ) » PR ( J » 3 1 ) » PR ( J » 46 ) t PR ( J > 62 ) > PR ( J » 77 ) 


1 1 PR ( J » 92 ) f PR ( J » 1 07 ) » PR (J»123)» PR (J»1Z8>» PR ( J » 154) »PR( J» 169) 
116 FORMAT (1211) 

103 FORMAT ( 12F6. 4) 

N=0 


DO 2 1=1.6 
DO 3 K=N+1 »N+15 
A( J.K ,, =A( J.N+1 ) 
C ( J . K ) -C ( J » N+ 1 ) 
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0040 

0041 

0042 

0043 
004* 
0045 
CO** 

0047 

0048 

0049 

0050 

0051 


0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 


0067 

0068 
0069 


PR<J*K)=PR<J»N+1> 

3 CONTINUE 

DO 4 K=N+ 1 6 » N+NDAY < I > 
A( J»K)=A< J»N+16) 

C< J»K)=C< J.N+16) 

PR< J»K>=PR< J.N+16) 

4 CONTINUE 

2 N=N+NDAY ( I ) 

40 CONTINUE 

READ( 1 1 107 > X« Y* RAi RB 
107 FORMAT < 4F10.5) 

CALL CLOSE < 1 > 
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c 

c 

c 



READ CLIMATILOGICAL DATA 
CARD GROUP It 

YEAR BEING MODELED *4 COLS 1-4 
CARD GROUP 2J 

DAILY DEGREE DAY <T> 16<12F5.2> COLS 1-72 

16 CARDS FOR EACH ZONE 
CARD GROUP 31 

DAILY PRECIPITATION 16<12F5.2> COLS 1-72 
16 CARDS FOR EACH ZONE 
CARD GROUP At 

DAILY SNOW COVER 16<12F5.4> COLS 1-72 

16 CARDS FOR EACH ZONE 
CARD GROUP 5* 

STREAM FLOW MARCH 31 F5.0 COLS 1-1G 

CARD GROUP 6J 

MEASURED STREAMFLOW 16<12F5.0) COLS 1-72 


CALL ASSI6N< 2 1 'CL I MAT * DAT ' 1 10. ' RDO ' ) 
READ <2» 110) YEAR- 
110 FORMAT < 14 ) 

DO 6 J= 1 » NZ 

6 READ <2 » 104 > < T < J » K ) » K=1 » 183 ) 

104 FORMAT ( 12F5. 2) 

DO 7 J=1»NZ 

7 READ<2» 104) <P< J.K) »K=1 » 183) 

DO 8 J=1 »NZ 

8 READ <2* 105 ) < S< J» K ) » K=1 » 183) 

105 FORMAT <12F5» 4) 

READ <2» 106) GN 
READ ( 2 * 106 ) ACTFLO 

106 FORMAT < 12F5.0) 

CALL I DATE < DATE < 1 ) . DATE < 2 ) » DATE <3 ) ) 

C 

WRITE ALL INPUT DATA 


DO 119 J= 1 » NZ 

WRITE<6» 121) NAME, YEAR, DATE. Jt < <M0( KN» K)»K=1»3)» KK=4 » 9 ) 

121 r ORMAT < 1H1 » //» 5X» ' RANGO/MARTINEC MODEL VERSION RCI-1.4 RUN FOR' 
120A4»//»20X, 'THE RUN IS FOR THE YEAR OF ' »I4» 

1' RUN OF MODEL MADE ' » I2» '/' , 12, '/' f I2» 

2//.30X, 'DAILY CLIMATOLOGICAL DATA FOR Z0NE',I2» 
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0070 

0071 

0072 

0073 

0074 

0075 

0076 


0077 

0078 
007? 

0080 

0081 

0082 

0083 


0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 


C 

C 

C 

C 

C 

C 

C 

C 


3// » 1 X » ' DAY ' »6<4X»3A4>4X» >* 
4/.4X»6< / DEGREE PREC. SNOW ')» 
5/*4X»6< * DAYS IN. COVER ')*/> 
DO 118 K=1 » 30 
K2=K+30 
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K3-K461 

K4»K+?1 

K5=K+122 

K6=K+153 

WRITE (6. 122>KtT( J.K) ,P< J,K)» S( J»K)»T< J»K2)»P< J.K2) tS< J»K2) 


1»T(J»K3) »P< J»K3) »S( J»K3) 

1 »T ( Jf K4> »P( J»K4 > » S< J»K4> » T ( J»K5) t P< J»K5 ) » S( J»K5) 

1»T( J.K6) tP(Jt K6 ) » S ( J » K6 ) 

122 FORMAT ( 1X» I2»2X»6( 1X.F5.2. 1X»F5.2» IX. F5.3.2X) ) 

118 CONTINUE 

WRITE <6* 126) T<J»61)»P<J»61)»S<J»6i)»T(J» 122) » P( J* 122) » S( J» 122) 
1»T( J»153)»P< J.153) ,S< J.153) 

126 FORMAT (1X»'31'» 22X» 1X»F5.2»1X*F5.2» 1 X»F5. 3» 2X» 20X 
1 »2<1X»F5.2» IX »F5 . 2t 1X» F5 . 3» 2X> ) 

119 CONTINUE 

WRI TE < 6 » 1 25 ) NAME » YEAR » DATE » < ( MO < KK » K ) t K= 1 » 3 ) » KK=4 » 9 ) » NDAY 
125 FORMAT ( 1H1 » //» 5X» 'RANGO/MARTINEC MODEL VERSION RCI-1.4 RUN FOR'. 
120A4»//»20Xt 'THE RUN IS FOR THE YEAR OF '.I4» 

1' RUN OF MODEL MADE 7 » I2» '/' » I2» '/' » I2» 

2// »30X» 'WATERSHED MELT AND RUNOFF COEFFICIENTS' t 
3// » 20X > 6 ( 3X » 3 A4 * 3X ) * 

4/»20X»6< ' 1-15 16 - ' t 12t IX) > 


DO 31 J=1 » NZ 
WRITE (6t 120) J 

120 F0RMAT<1H0»3X» 'DATA FOR ZONE 'tllf'?') 

WRITE < 6 » 123 ) A ( J » l)*A<J»16)rA(J»31) > A < J» 46) , A( J*62) , A< J» 77) » A( J»92 
1»A( J» 107) t A< Jt 123) »A< J»138) * A< J» 154) t A< J»169> 

123 FORMAT <6X> 'ME* T FACTORS '» 2X » 6< 3X, F5 . 3» 3X t F5 . 3 » 2X) ) 

WRITE <6»124)C(J»l)»C(J»16)»C(J»3l)»C<Jf 46 ) ,C< J.62) .C< J» 77) >C< J»92 
1 fC( J> 107) »C< J» 123) »C( J» 138) »C( J*154) »C< J»169) 

124 FORMAT <6X» 'RUNOFF COEF . ' » 2X» 6< 3X» F5. 3. 3X> F5 . 3 » 2X) > 

WRITE <6* 127 ) PR ( J* 1 ) » F'R ( J? 16) , PR < J» 31 ) , FR< J* 46) » PR < J> 62) » FR < J» 77) 

1 * F‘R C J » 92 ) » PR < J » 1 07 ) » PRC J» 123) t PR< J» 138) »FR< J» 154) »PR< J» 169) 

127 FORMAT <6X. 'PREC. METHOD' »2X»6<4X»Il»7XtIl»5X)) 

31 CONTINUE 

WRITE <6 » 128) 

128 FORMAT ( 1 HO » //// » 6X » ' NOTE I PREC. METHOD » ' » 

1/.20X»'0 = PREC. RUNOFF COMPUTED FOR NON-SNOWCOVERED AREAS ONLY'* 
l/» 20X » ' 1 * PREC. RUNOFF COMPUTED FOR ALL AREAS') 

ALL DATA HAS BEEN ENTERED 


BEGIN COMPUTATIONS 

DETERMINE IF PRECIPITATION IS RAIN OR SNOW 
AND IF RUNOFF OCCURS FROM PRECIPITATION OR- 
IS ACCUMULATED AS SNOW 


0096 

0097 


DO 10 J=*1,NZ 
SNOW * 0.0 
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0098 

0099 

0100 

0101 
0103 

0105 

0106 
0108 
0110 

0112 

0113 

0114 

0115 

0117 

0119 

0120 
0122 

0123 

0124 

0125 

0126 

0127 

0128 


0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 
0142 

0144 

0145 

0146 

0147 

0148 


DO 11 K-l » 183 
5C0=0.0 

sci=o.o ORIGINAL PAGE IS 

IF<F‘R( JtK) «£G.O) SCO=S(JpK) qc POOR OUAUTY 

IF<F‘R< J»K) . EQ. 1 > SC1«S<JpK> wr rvrw 4- 

R - T < J»K)*A< J»K) 

C CHECK FOR RAIN IN ZONE J 

IF<P< J» K ) . EQ. 0. ) GO TO 12 
C...... CHECK FOR SNOW IN ZONE J 

IF<T< JtK) .GT.2.0) GO TO 12 

C CHECK IF SNOWCOVER IS GREATER THAN 90 PERCENT 

IF<S< JtK) .6E.0.9) 60 TO 14 
C ACCUHULATE SNOW IN ZONE J 


SNOW = SNOW + P ( J » K > 
14 P< J? K)=0«0 
GO TO 11 


C IS THERE PRECIPITATION TO BE MELTED 

12 IF ( SNOW . EQ . 0 . ) GO TO 11 

C CHECK TEMPERATURE AND CACULATE SNOW TO BE MELTED 


IF<T<J»K).LE.O.) GO TO 11 
EP = R*<1.-SC1> 

IF<R. GT. SNOW) GO TO 13 
P<J»K)=P( J»K)+EP 
SNOW=SNOW-R 


GO TO 11 

C CACULATED RUNOFF IS GREATER THAN SNOW BEING MELTED 

13 P< J» K)=P< J * K ) + ( SNOW* < 1 . -SCI ) ) 

SNOW=0.0 

11 P< Jt K) = P<JpK)#(1.- SCO ) 

10 CONTINUE 
C 

C PRECIPITATION WHICH WILL RUNOFF HAS BEEN 

C CALCULATED FOR EACH ZONE 

C...... 

C BEGIN CACULATION OF DAILY STREAM FLOWS 

C 

IK) 15 K=1 » 183 
SUMRO =0.0 


SC = 0.0 
SA = 0.0 
DO 16 J=1 » NZ 

SC = SC + S< JtK)*AREA< J) 

SA = SA + AREA< J> 

RUNOFF < Jt K )=C < J» K ) #AREA< J ) *26 • 8889D0* <A(J»K)*T<J»K) #S( J t K ) +F'< J t K) 
16 SUMRO=SUMRO+RUNOFF< JpK) 

DR = RA + <RB*SC/SA) 

RECESS =X* < QN**Y > 

IF<K.EQ.l) F‘RUNOF=QN 
IF<K.E0.1) DRP = DR 

STRFLO <K) = < < 1 -DRP ) *PRUNOF+DR*SUMRO ) * ( 1 . 0-RECESS ) +GN*RECESS 

PRUNOF =SUMRO 

QN=STRFLO<K) 

DRP = DR 
15 CONTINUE 
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* 


014? 

0150 

0151 


0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 


0166 

0167 

0168 
016? 

0170 

0171 

0172 


0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 
0183 


..... WRITE COMPUTED RUNOFF FOR EACH ZONE 

DO 1? J=1 pNZ 

WRITE <6» 113) NAME » YEAR. DATE » J» ( <MO<KNt K) » K=1 * 3 ) » KK=4 .9 ) 

113 FORMAT (1H1.//.5X. 'RANGO/MARTINEC MODEL VERSION RCI-1.4 RUN FOR'# 
1 20A4 » // » 20X » ' THE RUN IS FOR THE YEAR OF '.14. 

1' RUN OF MODEL MADE ' t 12. '/' . 12. '/' t I2» 

2//. 30X. 'COMPUTED DAILY RUNOFF FOR ZONE'. 12. 

3//» IX » 'DAY . 6< 4X» 3A4 » 4X. ) t 
4/»4X»6< 1X» ' COMPUTED RUNOFF '.2X)» 

5/»4X»6<3X» ' CFS '*4X>»/> 

DO 18 K=1 » 30 
K2=K+30 
K3-K+61 
K4=K+91 
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K5=K+122 

K6-K+153 

WRITE(6* 1 14)K. RUNOFF< J.K ) » RUNOFF < J.K2) » RUNOFF ( J»N3) 

1 » RUNOFF < J » K4 > » RUNOFF ( J » K5 ) * RUNOFF < J » K6 ) 

114 FORMAT < 1X» I2»6<6XtF7. 1 »7X> > 

18 CONTINUE 

WRITE(6»115) RUNOFF ( J> 61 > * RUNOFF ( J. 122 ) »RUNQFF< J» 153 ) 

115 FORMAT <1X» ' 31 ' » 20X » 6X » F7 . 1 » 7X» 20X » 2< 6X » F7 . 1 » 7X> ) 

19 CONTINUE 

C WRITE COMPUTED STREAMFLOWS 

WRITE(6> 10?) NAME .YEAR. DATE. (<M0< J.K).K=lf3)f J=4»?) 

109 FORMAT ( 1H1 . // » 5X » ' RAN60/MARTINEC MODEL VERSION RCI-1.4 RUN FOR'. 
120A4.//.20X. 'THE RUN IS FOR THE YEAR OF '.14. 

1' RUN OF MODEL MADE ' . 12. '/ ' . 12. ' / ' . 12 . 

2///» 5X » 6 < 3X . 'STREAMFLOW F0R'3X>» 

3/» IX. 'DAY ',6(4X.3A4.4X. )» , 

4/.4X.6UX. 'COMPUTED' .3X» 'ACTUAL'2X>. 

5/.4X»6(3X» 'CFS' »7X. 'CFS' »4X) »/) 

DO 17 K=1 .30 

K2=K+30 

K3=K+61 

K4=K+91 

K5=K+1 22 

K6=K+153 

WRITE<6. 112)K.STRFL0<K) » ACTFLO<K) .STRFL0(K2) . ACTFL0(K2) . 
1STRFL0<K3) .ACTFL0(K3) »STRFL0(K4) . ACTFL0(K4) .STRFL0<K5) »ACTFL0(K5) 
2STRFL0 < K6 ) . ACTFLO < K6 ) 

112 FORMAT (IX .12. 6< 2X . F6 . 0 » 4X . F6 . 0 . 2X > ) 

17 CONTINUE 

WRITE(6. Ill) STRFLO< 61 ) .ACTFLO <61 ) »STRFLO< 122) . ACTFLO< 122) . 
1STRFL0< 153). ACTFLO <153) 

111 FORMAT < IX. '31 ' » 22X . F6 . 0. 4X » F6. 0. 22X » 2< 2X» F6 . 0 » 4X » F6 . 0. 2X) ) 

CALCULATE NASH-SUTCLIFFE GOODNESS OF FIT 

G=0.0 
H=0 . 0 

DO 200 1=1.183 
G=G+ACTFLO< I ) 

H=H+STRFLO< I ) 

200 CONTINUE 

GMEAN=G/1 83 . 
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0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 


0193 

0194 

0195 


0196 

0197 

0198 

0199 

0200 
0201 


E*0*0 

F=0»0 

DO 210 1-1 » 183 
E*E+ ( ACTFLO < I ) -GMEAN ) **2 
F«F+<ACTFLO<I)-STRFLO<I) >**2 
210 CONTINUE 

XNSR2= < E-F ) /E 

WRI TE ( 6 * 2 1 5 ) NAME » YEAR # DATE » XNSR2 

215 FORMAT <1H1*//*5X» ' R ANGO/MART I NEC MODEL VERSION PCX-1*4 RUN FOR'* 
i20A4»//*20X* 'SUMMARY OF THE MODEL RUN FOR THE YEAR OF '*14* 

1' RUN OF MODEL MADE ' » 12* * 12* '/' * 12#////* 

125X* 'NASH-SUTCLIFFE GOODNESS OF FIT METHOD - NSR**2 « '»F10.4»/> 

C* CALCULATE SEA1 ONAL DIFFERENCE 

PE=< <H-G)/H>*100.'. 

WRITE (6* 220 )G i H* 

220 FORMAT <25X' ACTUAL SEASON VOLUME = '#F12.3»' CFS-DAYS ' » //* 

*25X' COMPUTED SEASON VOLUME * '*F12.3*' CFS-DAYS'*//* 

*25X ' SEASONAL DIFFERENCE IN PERCENT - '»F10.2> 

C WRITE OUTPUT FILE FOR PLOTTING 

CALL ASSIGN<3» 'OUTPUT.DAT' » 10* 'NEW') 

WRITE <3* 30) STRFLO 
30 FORMAT < 1 2F6 • 0 > 

WRITE(3*30) ACTFLO 
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FORTRAN IV 


Storasae Map for Prosram Unit 5N0MEL 


Local 

Variables* .PSECT 

♦DATA* 

Size 

« 075274 

<15710. 

words) 


Name 

Type 

Offset 

Name 

Twee 

Offset 

Name 

Type 

Offset 

DR 

R*4 

075146 

DRP 

R*4 

075152 

E 

R*4 

075172 

EP 

R*4 

075132 

F 

R*4 

075176 

6 

R*4 

075156 

H 

R*4 

075162 

1 

1*2 

075062 

J 

1*2 

075056 

K 

1*2 

075064 

KK 

1*2 

075076 

K2 

1*2 

075100 

K3 

1*2 

075102 

K4 

1*2 

075104 

K5 

1*2 

075106 

K6 

1*2 

075110 

N 

1*2 

075060 

NZ 

1*2 

075054 

PE 

R*4 

075206 

PRUNOF 

R*8 

075000 

QMEAN 

R*4 

075166 

ON 

R*8 

075010 

R 

R*4 

075126 

RA 

R*4 

075066 

RB 

R#4 

075072 

RECESS 

R*8 

074770 

SA 

R*4 

075142 

SC 

R*4 

075136 

SCO 

R*4 

075116 

SCI 

R*4 

075122 

SNOU 

R*4 

075112 

SUMRO 

R*8 

075040 

X 

R*8 

075020 

XNSR2 

R#4 

075202 

Y 

R*8 

075030 

YEAR 

1*4 

075050 


Local and COMMON Arrays! 


Name 

Type 

Section 

Offset 

Si 

ze 

Dimensions 

A 

R*4 

Vec 

♦DATA 

025344 

007114 

< 

1830.) 

(5*183) 

ACTFLO 

R*4 


♦DATA 

063404 

001334 

< 

3 66.) 

<183) 

AREA 

R*4 


♦ DATA 

043574 

000024 

< 

10.) 

< 5) 

C 

R*4 

Vec 

♦DATA 

034460 

007114 

< 

1830. ) 

<5* 183) 

DATE 

1*4 


♦DATA 

074210 

000014 

< 

6.) 

< 3 ) 

MO 

1*4 

Vec 

♦ DATA 

074224 

000220 

< 

72.) 

< 12 * 3 > 

NAME 

1*4 


♦ DATA 

074070 

000120 

< 

40. ) 

< 20) 

NDAY 

1*2 


♦DATA 

074054 

000014 

< 

6.) 

< 6 ) 

P 

R*4 

Vec 

♦ DATA 

016230 

007114 

< 

1830. > 

(5*183) 

r R 

1*4 

Vec 

♦ DATA 

064740 

007114 

< 

1830.) 

(5*183) 

RUNOFF 

R*8 

Vec 

♦ DATA 

043620 

016230 

< 

3660. ) 

(5*183) 

S 

R*4 

Vec 

♦ DATA 

007114 

007114 

< 

1830. ) 

(5*183) 

STRFLO 

R*4 


♦ DATA 

062050 

001334 

< 

366. ) 

(183) 

T 

R*4 

Vec 

♦ DATA 

000000 

007114 

< 

1830. ) 

(5*183) 


Subroutines* Functions* Statement and Processor-Defined Functions! 


Name 

Type 

Name 

Type 

Name 

Type 

Name 

Type 

Name 

ASSIGN 

R*4 

CLOSE 

R*4 

I DATE 

1*2 
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appendix b 


Summary of Published Climatological Data 
for National Weather Service Stations at: 
Pagosa Springs 
Wolf Creek IE 
Del Norte 
Hermit 7ESE 


April - September, 1973-1980 


SUMMARY OR Cl. IMATULtlfi TCAL PATA FOR THE YEAR 1973 
FOR THE MONTH OF ATRIL 
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Adjusted Streamflow Records 

Rio Grande 
1973-1980 

South Fork 
1973-1980 
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MEASURED STREAM/ LOWS ADJUSTED TO ACCOUNT FOR RESERVOIR STORAGE 


40 

50 

ft 

80 

8S 

90 

100 

18 

130 

140 

ISO 

II 

190 

200 

210 

220 

230 


REM mmmmmtmiMNusTmmtmmmmtmttmtmt 

WtfSHA 

RESERVOIR STORAGE AMD RELEASES. 

mmmtmm initialize variables mttmmmtmm 


REM 

REM 


REM 

REM 


BIN RISK 

dim s<B): 

DIM W35): 
BIN V 3SK 
PIH HIT 
BIN T(! 

sink: 

BIN Z(2 
BIN A(35)J 

y=o: 

rim>: 

si=o: 

ci=o: 


REN (R.C. RES. STORAGE FOR 
(S.N. KS. STORAGE FOR 


BAY OF 
DAY OF 


MONTH) 

MONTH) 


REM 


350 

349 


REN 

REM (CONT, RES. STORAGE FOR EACH DAY OF MONTH) 

REN (OBSERVE! CAGE READING ATjgL NORTE) 

REN (CHANGE IN R.C. RES. STORAGE) 

REN (CHANGE IN S.N. RES, STORAGE) 

(CHANGE IN CONT. RES, STORAGE) 

(SUN OF ALL THREE RESERVOIR STORAGES 
CONVERTER TO ACRE-FEET/DAY) 

REN (ADJUSTER GAGE REARING AT BO. NORTE) 

REN (YEAR) 

REN (RIO GRANDE RESERVOIR STORAGE CRAY N411) 

•Td I (SANTA MARIA RESERVOIR STORMS [DAY Ntll) 

REN (CONTINENTAL RESERVOIR STORAGE [DAY NU1) 

KM (SUNJf RESERVOIR CHANGES FOR BAY 0) 
m, (tfJNlftff BAYS IN MONTH ) 

REM (BAY » MONTH SUBSCRIPT) 

REM (PRINT FORMAT CHARACTER) 

REM (PRINT FORMAT CHARACTER) 

REM (FORMFEED CHARACTER) 

REM (MONTH) 

REN (CORRECT DffliT dEOC) 

ISM (IMffiER OF BAYS WITH ZERO STORAGE [NO DATA!) 

REN (AVERAGE) 

m (SUBSCRIPT FOR FIRST BAY WITH 0 STORAGE [NO BATA]) 

aumimmtt set up input bata mtmmmut 

INPUT ‘HHAT MONTH* I M 

w* 

I»UT ‘CORRECT* PUSH RETURN* INCORRECT, PUSH -1 AND RETURN* II 
IF IOQ GOTO 380 

D0TITT * • 


A$=‘ZZZZZZZZZ M . 
B4=*ZZZ*; 
F$=CHARK 12)5 

w: 

i=o: 

k=o: 


490 REM 
700 

ft 

710 
720 
730 
740 


INPUT *UHAT YEAR* I Y 
1=0 

PRINT Y 

INPUT ’CORRECT* PUSH RETURN! INCORRECT, PUSH -1 AND RETURN* II 
IF IOC GOTO 440 
PRINT * * 

PRINT * * 

INPUT *HQU MANY DAYS IN MONTH* iN 
PRINT N 
1=0 

INPUT 'CORRECT, PUSH RETURN* INCORRECT, PUSH -1 AND RETURN’ II 
IF IOO GOTO 540 
PRINT* • 

PRINT* * 

PRINT *RIO CRANK RESERVOIR STORAGE FOR BAY ONE OF NEXT MONTH’ 
PRINT *IF NO BATA, ENTER VALUE FOR LAST BAYOF THTS MONTH* 

INPUT R1 
PRINT Rl 
1=0 

INPUT ‘CORRECT, PUSH RETURN! INCORRECT, PUSH -1 AND RETURN’ il 
IF IOO GOTO 420 
PRINT * ‘ 

PRINT ‘ ‘ 

PRINT ’SANTA MARIA RESERVOIR STORAGE. FOR BAY ONE OF NEXT MONTH’ 
JRjjjT *jF NO BATA, ENTER VALUE FOR LAST DAY OF THIS MONTH* 

PRINT SI 
1=0 

INPUT ‘CORRECT, PUSH RETURN! INCORRECT, PUSH -1 AND RETURN* !l 
IF IOO GOTO 700 
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750 


REM 


PRINT ' 
PRINT • 


781 

782 


$ B 


PRINT 'CONTINENTAL RESERVOIR STORAGE FOR DAT ONE OF NEXT MONTH* 
PRINT 'IF NO DATA. EN« Jt VALUE FOR LAST DAY OF THIS MONTH* 

INPUT Cl 
NT Cl 


611 

820 


REM 


930 

940 

8 

970 

980 

loot 

1916 

1020 

1030 

118 

1065 

1070 

1080 

1090 

HOC 

1110 

1130 

1140 

1150 

1155 

118 

1180 

1190 

m 

1215 

1221 


INPUT ‘CORRECT! PUSH xETURNI INCORRECT* PUSH -1 AND RETURN* 11 
IF JOO GOTO 780 
‘ NT * * 

NT * • 

'PRINT 'SUM OF CHANGES IN RESERVOIR STORAGE FOR DAT 0?* 

PRINT *(IF FIRST MONTH OF TEAR, ENTER 0)* 

INPUT Z1 
PRINT Z1 
M 

“"‘I. I RETURN} INCORRECT* PUSH -1 AND RETURN* II 

COT 

PRINT * * 

PRINT * * 


ns n»" 


REM 

REM tmitttt SET UP RESERVOIR STORAGE AN) OWNS IN STORAGE ARRAYS Mill 

rem timtmmmmu for rio grande reservoir ttmmmmmmm* 

REM 

FORM TON 
PRINT * * 

PRINT * * 

PRINT 'RIO CR/WE RESERVOIR STORAGE FOR*,** '}D4,*' '}Ji'» 'WY* 
PRjjjT 'jnjO DATA, ENTER 0* 

PRINT P45* , }J}*, , *T»'= , *R< J) 

1=0 

INPUT •CORRECT, PUSH RETURN} INCORRECT, PUSH -1 AND RETURN* fl 
IF IOO GOTO 1030 
NEXT J 
K=0 

FOR J=1 TO N 

IF WJH) GOTO 1230 
KIMKJMMKJ) 

IF KOO GOTO 1180 

¥U GOTO 1255 
A9=( R( HFK )-W M-l ) )/< Kf 1 ) 

FOR L=< H-l ) TO (NFN-l) 

H(U=A9 
NEXT L 
K=0 

GOTO 1160 


124 
125C 
1255 
1265 
IK- 
1285 
1295 
1300 
1310 
1320 
1330 
1340 
1350 
1365 
1370 


IF K=0 THEM H=J 
K=Ktt 
GOTO 1160 
WNttHKN) 


1400 

1410 

1430 

1440 

1450 

1455 

1460 

1470 

1480 

1490 

iSI 


tttmttt SET IT RESERVOIR STORAGE AMD CHANGE IN STORAGE ARRAY Wit 

ren mmmmmmmtt for santa maria reservoir tomtit mum 
REM 

FOR M TO N 
PRINT * 1 
PRINT 1 * 

PRINT 'SANTA MARIA RESERVOIR STORAGE FOR'}* *}&$}• •IJ**, *iYi'?* 
PRINT *IF m DATA, ENTER 0* 

INPUT S (J) 

PRINT Hi* '}J}*, 'JYJ^'iSU) 

1=0 

RETURN,* INCORRECT, PUSH -1 AND RETURN'} I 

NEXT J 
K=0 

FOR J=1 TO N 

IF SUM GOTO 1530 
TUWU4D-9U) 

IF KOO GOTO 1480 
NEXT J 

IF K=9 GOTO 1555 
AN S<ltt*)-SCIH ))/<*«> 

FOR L4H-IM0 <MK-1) 

T(L)=A9 
NEXT l 


.071 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2165 

90 

2170 

2)80 

2190 

2200 

2210 

2220 

2230 

2240 


C0?0 1460 

COTO 1460 
T(N)=S1-S(N) 



1850 

1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
‘flto 
w) 
050 
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ItltttU SET UP RESERVOIR STORAGE AW 0 ^ IN STORAGE ARRAY mill 
m i«mmmmm«FOR continental reservoir tmmmtmmm 
REN 

FORM TO N 
PRINT • ■ 

PRINT * * 

PRINT 'CONTINENTAL RESERVOIR STORAGE FOR*?* •IW* '»J?', 'TY?'?' 
PRINT 'IF NO DATA* ENTER 0* 

i win cm 

*»CU) 

INCORRECT, PUSH -1 AW RETURNS I 

. J WUIO 1MO 

NEXT J 
K=0 

FOR 2=1 TO N 

IF CU)=0 COTO 1830 
MJKMHHIJ) 

IF KOO COTO 1780 
NEXT 4 

IF K=0 COTO 1855 
A9=UXN«HXN-1))/(XU) 

FOR L=<N-1) TO (WK-1) 

Kl)=A9 
NEXT L 
K*0 

WTO 1760 
IF K=0 THEN H=J 
Mitt 
WTO 1760 
WN)=C1-C<N) 

mmmttm set up observes streanflow data mmmmm 

REN 

FOR J=1 TO N 
PRINT • • 

PRINT • • 

PRINT 'ENTER OBSERVED STREAfflOW DATA FOR'?* *?Bt,*' *?J?', 'JYJ'T* 
INPUT QU) 

PRINT D9»' "Ml'» •.‘Y»'='}Om 
1=0 

DWT 'CORRECT, PUSH RETURN? INCORRECT* PUSH -1 AND RETURN* ?I 
IF IOO WTO 1920 
IIEXT 4 
REN 
REN 

RENOTW COHPUTF SUN OF ALL THREE RESERVOIR STORAGE CHANGES (Z* AND till 
mum CONVERT TO C.F.S. PER DAY (DIVIDE A.F. DY 1.98347) tIUlimm 
REN 

FOR J=1 TO N 

Z(J)=<H(JWT(J)4MJ))/1. 98347 
fc(T J 
REN 

RENHWI ADJUST OBSEPVED KL NORTE READING (0) FOR THE SUN OF THEIittl 
RENMW RESERVOIR STORAGE CHANGES (ADD ONE DAY TRAVEL TIHEMttUt 
REN 

A( 1 )=0( 1 H71 
FOR J=2 TO N 
AWWNJHKM) 

€XT i 
REN 

ren mmmtmmm print table ttmmutmmmut 

REN TURN ON PRINTER 
P^INT CHARK17) 

REN 

REN FORNFEEB 
PRINT FI 
REN 

PRINT • 

PRINT • 

PRINT ' 

PRINT • 


ADJUSTED STREMF1.W AT DEL NORTE, COLORADO! RTO GRANDE RJ9ER' 
(ADJUSTED FOR RIO CRANDE, SANTA N»IA» AND CONTINENTAL RESERVOIRS)' 


9 


IMI' »SY 


2250 PR 
2260 PR 
2270 PR 
2280 PR 
2290 PR 
230# PR 

« 


n ■ 

|i 


INT 

INT • BATE 

IMT * 

IHT * ■ 

2330 FORJ*l TON 
“ IRT * • 

IRT mTtJ.H)}* 


RIO CRANDE 
STORACECHAMCE 


SANTA MARTA 


CONTINENTAL 


CAGE READING ( BEL WV 


STORAGE CHANCE STORAGE CHANCE JUH CHANCE < CFS/BAY ) W3EWEI AFUySTEB* ( 


2340 pr: 

2350 PRL 
Z(J).A9)f 

2$0 ^PRinf F4 
2110 PRINT CHARM 19) 
REACT 


FHT( « 4 )»Af )i RfT( m ). At )i FHT< S* J )» At M* • I FMT< T< 4 ).A4 )} FHT( C' 4 )» At >• 1 ' • i FW K ■» )»A4 )* ' 


i FMT(OU)iAI)}* '{ FHTtm J)fA$ ) 
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APPENDIX D 

Snqv-Cover Depletion Curves 
Rio Grande River above Del Norte 

1973-1981 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1970 
FOR THE MONTH OF SEPTEMBER 
DEG. DAY - < TMAX+TMIN >/2. 0 - 32. 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1977 
FOR THE MONTH OF MAY 
DEO. DAY - < TMAX+TMIN)/2 • 0 - 32. 
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DEO. DAY - (TMAX+TM2N)/2.0 - 32. 


ORIGINAL PA‘3£ StS 
OF POOR QUALITY 



N 

9 


rv 

N 


n 

o 


n 

K> 


J N 
h- IV 

a 


A SUMMARY OF MODEL DATA FOR THE YEAR 1977 
FOR THE MONTH OF JULY 
DEO. DAY - <TMAX4TMIN>/2.0 - 32. 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1977 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1*78 
FOR THE MONTH OF APRIL 
DE8. DAY - <THAX+THIN)/2.0 - 32. 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1978 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1*78 
FOR THE MONTH OF JUNE 
DEO. DRY • <TMAX+TNIN>/2»0 - 32. 
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SUMMARY OF MODEL DATA FOR THE YEAR 1978 
FOR THE MONTH OF AUGUST 
8, DAY - <TMAX+THIN>/2.0 - 32. 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1*7* 
FOR THE MONTH OF JUNE 
DEB. DAY - (TMAX+THIN)/2. 0 - 32. 
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A SUMMARY OF MODEL DATA FOR THE YEAR 1979 
-FOR THE MONTH OF AUGUST 
DES. DAY - (TMAX+THIN)/2.0 - 32. 
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A SUMMARY OF MGHEL DATA FOR THE YEAR 19*0 

FUR- THE— MOUTH -QL AFfcU 

DEG, DAY » < TMAX + TMIN > / 2 , 0 - 32 . 
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A SUMMARY OF MODEL DATA FOR THE YEAR 19*0 
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u> 

’ w 
* X 

u 

ui 

.t 


N •» » 

; ao 


U a 
u 
* 

(9 

U 

a 


<x 

a 







A SUMMARY OF MODEL DATA FOR THE YEAR 1960 
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DEG. DAY - < TMAX+TMIN ) /2 . 0 - 32. 
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A SUMMARY OP MODEL DATA TOR THE YEAR I960 

FOR— T HE— MON TH—OR JULY- 

DEG. DAY » < TMAX+TMIN)/2 .0 - 32. 
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SUMMARY OF MOi'EL DATA FOR THE YEAR 1980 
FOIL THE— MONTH OF AUGUST 
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